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Abstract 

We show that the gravitational binding energy of a neutron star of a given mass is correlated with 
the slope of the nuclear symmetry energy at 1-2 times nuclear saturation density for equations of 
state without significant softening (i.e., those that predict maximum masses Mmax > 1.44Mq in line 
with the largest accurately measured neutron star mass). Applying recent laboratory constraints 
on the slope of the symmetry energy to this correlation we extract a constraint on the baryon mass 
of the lower mass member of the double pulsar binary system, PSR J0737-3039B. We compare 
with independent constraints derived from modeling the progenitor star of J0737-3039B up to and 
through its collapse under the assumption that it formed in an electron capture supernova. The 
two sets of constraints are consistent only if L < 70 MeV. 

PACS numbers: 21.65.-f,21.65.Ef,21.65.Mn,26.60.Kp,97.66.Gb,97.60.Jd,98.80.-d 
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I. INTRODUCTION 

The double pulsar system PSR J0737-3039A/B consists of two radio pulsars with 

masses and spin periods 1.3381±0.0007Mo, 22.7ms (pulsar A) and 1.2489±O.OOO7M0, 2.77s 
(pulsar B) in a compact 2.4h orbit with a relatively low eccentricity of e = 0.088. Precision 
observations of the system are providing excellent tests of a wide range of fundamental 
physics, including stringent tests of gravity in the strong field regime [3|]. J0737-3039A/B 
also provides insight into the evolutionary paths of binary systems that produce double 
neutron star systems (DNSs). 

DNSs originate in binaries where, necessarily, both stars are initially massive enough to 
undergo supernova (SN) explosions. We shall refer to the supernovae that created pulsars A 
and B as SNA and SNB respectively. It has been shown 4|-[6|] that the possible evolutionary 
paths of the binary lead to the immediate precursor of the DNS being a helium star (HeS) 
- neutron star (NS) system. 

The low mass of pulsar B and low eccentricity of the orbit (e = 0.088; estimated initial 
eccentricity eg ~ 0.11-0.12 [1]) led to the suggestion that it was formed, not by the collapse 
of an ton core, bnt by the collapse of an ONeMg core initiated by electron captures (e- 
captures) onto Mg; that is, an e-capture supernova [q. The gravitational mass of pulsar 
B is Mq ^ 1.25Mq, and its gravitational binding energy may be roughly estimated to 
be ^ OmAiMc/MQfMQ = O.UMq leading to an estimate of the baryon mass of the 
neutron star of Mb ~ 1.38M0. This is close to the estimated mass at which an ONeMg 
core becomes unstable to collapse . Such a supernova would be relatively symmetric as 
the timescale for the explosion is expected to be significantly shorter than the timescale for 
large instabilities to develop jojjlll: this would explain the low eccentricity of the system. 
Assuming a symmetric SN, the eccentricity implies that the mass expelled from the HeS 
progenitor during SNB was AM ^ eo(MA + Mb) = O.32M0 and the corresponding total 
mass for that HeS was ~ l.SGM©. Such a low mass helium star would form in a HeS-NS 
binary in which the initial HeS mass was ~ 3M0 {4,5, placing it in the mass range of helium 
stars expected to end their life in an e-capture SN [7|]. 

Additional evidence for a relatively symmetric supernova comes from the low value mea- 
sured for the transverse velocity of J0737-3039A/B, ^ 10 km s~^, which makes it statistically 
unlikely that SNB provided a large kick to the system 121]. This conclusion is additionally 
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supported by the system's position close to the galactic plane [l3|, the low alignment angle 



between the spin of pulsar A and the orbital angular momentum of the system 12|, and 
the stability of pulsar A's pulse shape profile 14 1. 

Under the hypothesis that pulsar B formed from an e-capture SN, a constraint on the 
neutron star equation of state (EoS) can be derived from modeling the progenitor's evolution 
up to and through the supernova jo]. As mentioned above, the progenitor ONeMg core 
becomes unstable to e-capture and collapses at a certain mass. Podsiadlowski et al jo] 
(hereafter referred to as article I) estimated that mass to be in the range 1.366 Mq < M < 
1.375M0. Under the assumption that negligible mass (< IO^^Mq) is lost from the core 
during the SN explosion, this then translates into the range for the baryon mass of pulsar 
B Mb. The gravitational mass is measured to be Mq = 1.2489±0.0007 Mq] the difference 
in the two masses gives the gravitational binding energy in units of mass BE = Mb — Mq. 
Note that since the gravitational mass is fixed at the measured mass of pulsar B, the binding 
energy and baryon mass of the neutron star are equivalent quantities. A given EoS predicts 
a certain binding energy for a neutron star of a given mass; therefore the constraints on Mb 
from modeling and the observed Mq together constrain the EoS under the e-capture SN 
hypothesis. Of the EoSs tested in article I, those based on field theoretical descriptions of 
nucleon-nucleon interactions appeared to give too small a binding energy to be consistent 
with the constraint. A later study using different field theoretical models found the same 
systematic under-binding ISj; however it is too early to say whether this is an intrinsic 
property of EoSs of that type. EoSs with different theoretical underpinnings showed no 
systematic under- or over-binding with respect to the original constraint, giving instead a 
scatter of binding energies 1.33M0 < Mb < 1.38M0. A later simulation of an e-capture SN 
from a 1.377M0 ONeMg core results in a neutron star of baryon mass 1.36 ± O.OO2M0 j^, 
indicating core mass loss of up to 1O~^M0. Such a mass loss shifts the constraint to lower 
binding energy, making it consistent with some of the field-theoretical models of the EoS. 

A pertinent question to ask is: can we identify specific properties of the neutron star 
EoS that correlate with the binding energy BE of a neutron star of a given gravitational 
mass? That is, what EoS property determines whether it will fall within the constraints 
imposed by modeling the e-capture scenario? In this article we shall demonstrate that such 
a property exists for EoSs of neutron stars that predict a maximum mass above 1.44M0. 
That property is the first derivative of the nuclear symmetry energy with respect to density 
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at nuclear saturation density no, commonly denoted L. Furthermore, given the correlation 
between L and BE, constraints on L from a variety of terrestrial nuclear experiments allows 
us to place an independent constraint on the binding energy of pulsar B, and hence ask: is 
the modeling of the baryon mass of pulsar B under the assumption of e-capture SN consistent 
with nuclear laboratory data? 

The article is organised as follows: in section [TTl we introduce the nuclear symmetry 
energy and its relation to the neutron star EoS as well as the current constraints on its 
density dependence from nuclear experiments. In section IIIII we summarize the wide range 
of EoSs that we will use to test for a correlation between symmetry energy and baryon mass, 
and in section IIVI we demonstrate the correlation between L and the neutron star binding 
energy, and use this correlation together with the experimental data to constrain the binding 
energy of pulsar B. In section |V] we discuss the implications of this result and in section IVl] 
conclude. 



II. THE NUCLEAR SYMMETRY ENERGY 



For completeness, we recall the relation of the nuclear symmetry energy to the neutron 
star equation of state. 

The neutron star equation of state, in the region of neutron, proton, electron, and muon 
(npe/i) matter is determined by the equation of state of infinite nuclear matter, which may 
be written as an expansion of the energy per nucleon in terms of the isospin asymmetry 
parameter a = 1 — 2yp, where yp = rip/n is the proton fraction and rip, n the number 
density of protons and the total baryon number density, respectively: 

E{n, a) = -EsNM(n) + a'^Esym{n) + ... (1) 

-E'snm(^) is the energy per nucleon of symmetric nuclear matter (proton fraction = 1/2 
or a = 0) and the nuclear symmetry energy Esymin) is defined 



_ , d^Ein,a) 



(2) 

«=o 



The expansion is in even powers of a only because of the approximate isospin symmetry 
of the strong force. a'^Esym{n) is the energy required to convert a certain fraction a of 
protons into neutrons from symmetric nuclear matter at a given density n. The parabolic 
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approximation, in which Eq.(IT]) is truncated at second order, is often used, and is generally 
accurate to high isospin asymmetries (a ~ 1) (although some properties of neutron star 
matter such as the crust-core transition density can be sensitive to higher order corrections 



16[); in that case, from Eq.([T]) the symmetry energy may be written Esy^{n) = i?pNM— -E'snm 
where -Epnm = E{n,a = 1) is the energy per nucleon of pure neutron matter. 

A fundamental property of symmetric nuclear matter is its saturation at a density of 
0.14 ^ riQ < 0.17fm~^ (equivalent to a mass density of ~ 2.5 x 10^^ g cm~'^), and an energy 
per particle of Eq = Es^^uino) ~ — 16MeV. It is customary to expand both the energy of 
symmetric nuclear matter and the symmetry energy in density around no. Defining the 
parameter x to be x = (n — no)/3no, those expansions can be written 

Estiuin) = EsNuino) + \Kox^ , (3) 

-Esym(n) = Esym(no) + LX+ ^KsyuiX"^ , (4) 

where the parameter L characterizes the slope of the symmetry energy at nuclear saturation 
density 



dE 



dE, 
io- 



3m ^^^^"^ 



dx ^ " n„ ' 

and Kq and -ft'sym characterize the curvature of the symmetric nuclear matter energy and the 
symmetry energy respectively. Note that, because symmetric nuclear matter has a minimum 
in energy at nuclear saturation density (by definition), there is no first order term in Eq.([3]). 
The pressure of infinite nuclear matter is 

PKa) = .^^. (6) 

Thus, at saturation density, the pressure of nuclear matter can be written 

P{nQ,a) = ^a^noL, (7) 

within the parabolic approximation. Therefore, at densities around saturation, the dominant 
contribution to the pressure of nuclear matter comes from the symmetry energy. 

Constraints on the magnitude and density dependence of the nuclear symmetry energy 
have been obtained in analyses of a variety of nuclear experimental data [l^, Q. Recent 



modeling of isospin diffusion in heavy ion collisions involving Sn and Sn extract con 



straints of 62< L <107 MeV 
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20| and 45< L <103 MeV 18| respectively. The extraction 



of the range of L is model dependent; the two ranges quoted here, whereas they overlap 



significantly, are from two different transport model ana^ 



from multifragmentation reactions yields L ~ 66MeV 2l|. A constraint from the analysis 



yses. A study of isoscaling data 



of pygmy dipole resonance (PDR) data gives 27< L <60 MeV 22| and from the analysis of 
;he surface symmetry energies of nuclei over a wide range of masses gives 75 < L <115 MeV 



231] . Finally, the measurement of neutron skins of a wide mass range of nuclei has led to an 



estimate of 25< L <100 MeV |24|. 

A conservative constraint, therefore, would be the full range given by the above studies 
25< L <115 MeV, which we shall call LI. As an illustrative example, we shall also examine 
the implications of the constraint of 62< L <107 MeV [lol, [2^ which we shall call L2. 



III. NEUTRON STAR EQUATIONS OF STATE 



A given nuclear matter equation of state completely specifies the EoS of uniform npe/j 
neutron star matter (the leptonic contribution to the EoS can be expressed in terms of 



the symmetry energy through the /3-equilibrium condition [17|, l25|). In Sec. |TT]we defined 
a number of parameters that characterize the magnitude and density dependence of the 
energy of nuclear and, hence, neutron star, matter: _Eo, -Esyml'^-o), -^5 -K^o, -^sym, etc. The 
sample of EoSs we use in this article span a representative range of these parameters. Our 
one constraint is that the EoSs used must predict a maximum neutron star gravitational 
mass of Mg > IAAMq in line with the maximum accurately measured neutron star mass at 
the current time {26 ]. 

An EoS is constructed using an underlying model of the nucleon-nucleon interaction and 
a given procedure to calculate the interactions of neutrons and protons in medium. We shall 
use EoSs with the following theoretical underpinnings (because we are using a large number 
of EoSs, we give mainly references to works that studied certain subset of them; individual 
references can be found in those works). 
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A. Phenomenological 



We use EoSs derived from 63 parametrizations of the nonrelativistic phenomenological 
Skyrme-Hartree-Fock model. They consist of those classified as group 1 and group II in 
that satisfy our NS mass requirement as well as the BSk6-16 parametrizations j28|. Also 
included are three phenomenological models used in article I: BPAL21, BPAL31 and BAL. 

B. RMF 

We use 29 EoSs based on a relativistic mean-field description of the nuclear interaction. 



We include all 22 used in a previous study [29| and we direct the reader there for details 



of those interactions. They include non-linear models, density-dependent meson-nucleon 



coupling models and point coupling models. We also include the KVR, KVOR 30| and D'^C 



31| interactions. Finally, we use four EoSs from article I: GLEN210,GLEN240,GLEN300 
and GLENHYB, the first three of which include hyperons and the fourth quark matter in 
the NS core. The latter four we shall collectively refer to as RMF-Hybrid models. 

C. Microscopic 

Six EoSs derived from potential models (based on realistic nucleon-nucleon potentials 
fit to scattering data together with a phenomenological three-body force) are used. BJ, 
FPS,APR,WFF and MOSZ [6] are calculated using non-relativistic methods whereas the 
DBHF EoS [32] is derived using a relativistic method. BJ, APR, MOSZ and WFF include 
hyperons at higher densities. 

D. MDI 



Also a phenomenological model, the modified Gogny interaction MDI [33| is used. We 
distinguish this from the other phenomenological EoSs because it was used in the extraction 
of the constraint on L, L2, from the isospin diffusion phenomenon 
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20| and because 



MDI contains a continuously variable parameter x that directly relates to the slope of the 
symmetry energy L. Each value of x gives a new symmetry energy dependence, and con- 
sequently a different EoS, while leaving the symmetric matter EoS the same. This allows 
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us to investigate the effect of varying L while keeping all other EoS parameters fixed at 
particular values. We generate EoSs by varying x between +0.6 and —2.0 (L = 32MeV and 
L = 153MeV respectively). 



E. QMC 



We use the two EoSs QMC700 and QMC7r4 calculated from the quark-meson coupling 
model 34], a relativistic theory in which nucleons are treated as bags of quarks that couple 
to meson fields. Hyperons are included self-consistently in these two EoSs. 



Neutron star models are calculated by solving the general relativistic equations of hydro- 
static equilibrium (the TOV equations); pulsar B is rotating slowly enough that corrections 
because of rotation are negligible. We take the gravitational mass to be fixed at 1.2489M0. 
The baryon mass is simply the number of baryons that constitute the star, A, multiplied 
by the mass per baryon ttlb, Mb = ArriB- An important point to note is that to equate 
that baryon mass with that of the progenitor core in the e-capture scenario, the mass per 
baryon must be taken to be me ~ 931.5MeV/c^, the value appropriate for a ONeMg core, 
rather than the nucleon mass ~ 939 MeV/c^. Taking the latter would give an ~ 0.8% shift 
in our baryon mass relative to that obtained by modeling the progenitor collapse, which 
is a significant difference in our context (e.g., the estimated uncertainty in baryon mass in 
article I is, O.OIIMq, or 0.8% relative to the total baryon mass.) 



IV. DEPENDENCE OF THE GRAVITATIONAL BINDING ENERGY ON THE 
SLOPE OF THE NUCLEAR SYMMETRY ENERGY 

It has been demonstrated that there exists a scaling relation between the radius of a 
neutron star R and its internal pressure P{n) sampled at a baryon number density n = 1-3 
hq. Empirically, it was found that for a wide range of neutron star EoSs which predict a 
maximum neutron star mass > 1.55Mq the relation is RP{n)~" ^ const., where a ~ 0.25 
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35| . The constant is dependent on the density at which the pressure is sampled, and is 
weakly dependent on the mass of the star. It was shown that the correlation is stronger the 
higher the density at which P was evaluated and the lower the mass of the star. We shall 
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evaluate the pressure at saturation density no- For the EoSs considered in this work, we 
plot i?P(no)~°'^^ versus R in Fig. [T| we find RP(n)~^'^^ ~ 9.5 ± 1. The scatter is somewhat 



larger than that found in the previous studies 25|, |35| , mainly because we allow EoSs with a 



lower maximum mass (1.44Mq rather than I.SSMq) and evaluate the pressure at a relatively 
low density (although these factors are compensated to a small extent by the relatively low 
mass of the neutron star). 

The above correlation can be expressed P{no) oc R^. As discussed in Sec.[Tll the dominant 
contribution to the pressure of neutron star matter at densities around saturation is from 
the derivative of the nuclear symmetry energy L, as demonstrated in Fig. [2] for our EoS 
selection. Indeed, it was pointed out that the R — P correlation translates into an R — L 



correlation 



35| : L oc -R^. Let us define the compactness parameter (3 = GMq/Rc^. At a 



constant gravitational mass, we can thus express the correlation L oc 

The binding energy of a neutron star relative to its gravitational mass, BE/Mq, is a 
function of compactness (3 (e.g., in Newtonian gravity it is simply BE/Mq = 3GMq/5R = 
0.6/3) where BE is in units of mass. In Fig. [3] we plot BE/Mq against /3 for our collection 
of EoSs. The analytic Newtonian relation is also plotted, and gives a good fit to the binding 
energy of the relatively low mass neutron star. 

To summarize: we have demonstrated for our model neutron star of gravitational mass 
Mq = 1.2489Mq, and for a wide range of EoSs, that the binding energy correlates roughly 
linearly with the compactness of the neutron star /3 which in turn correlates strongly with 
the slope of the nuclear symmetry energy L predicted by the EoSs. Therefore we expect to 
see a correlation between L and the binding energy of the neutron star, and consequently the 
baryon mass given that the gravitational mass is fixed. In Fig. H] we plot L versus the baryon 
mass of our model neutron star for all EoSs considered; such a correlation is immediately 
apparent. 



V. DISCUSSION 

Fig. m shows that, for a wide range of EoSs with widely varying density dependences 
and calculated using a number of different theoretical models of the nuclear interaction, 
a region in L — Mb space is preferentially picked out. We shall examine the implications 
of this correlation in light of the two e-capture constraints: Mb in the range 1.358-1.362 
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Mq [9| and Mb in the range 1.366-1.375 Mq [6!], which we shall refer to as Ml and M2 
respectively. In Fig. |U the constraints from nuclear experiment LI and L2 are represented 
by the dark and light shaded horizontal regions, respectively, and the two constraints Ml 
and M2 are defined by the dashed and dash-dotted vertical lines, respectively. For the two 
sets of constraints to show consistency we demand that the region of L — Mb space over 
which they intersect coincide with part of the region populated with EoSs. 

Let us first examine how the constraints on L, LI and L2, constrain Mb. The most 
conservative range of L, LI, is consistent with both Ml and M2; the constraint L2 coming 
solely from the analysis of isospin diffusion in heavy ion collisions with the MDI interaction 
is not consistent with M2 and only marginally consistent with Ml. One can see that values 
of L below ~ 60 MeV are consistent with a broad range of Mb, from ~ 1.35Mq up to 
at least ~ I.4OM0; Mb is only weakly dependent on L in this range. A value of L < 60 
MeV would therefore be consistent with the e-capture scenario, but not provide particularly 
stringent constraints on Mb while other properties of the NS EoS remain unconstrained. A 
value of L above 60 MeV would constrain Mb to within ~ O.O2M0 at values smaller than 
the constraint M2 and only marginally consistent with Ml. 

One can also consider how the current theoretical constraints from the e-capture SN 
scenario. Ml and M2, constraint L. Taking Ml and M2 together (giving a conservative 
range of Mb of 1.358 - I.375M0), the range of L is constrained to be L < 70MeV. The 
e-capture scenario thus appears consistent with the lower range of L values given by current 
experimental constraints. Looking at those constraints individually, only the range extracted 
from the analysis of the surface symmetry energies of nuclei 75 < L <115 MeV 23|] is 
completely inconsistent with the e-capture scenario. 

We finish our discussion with some caveats and remarks. Firstly, we are assuming that 
pulsar B is a neutron star; that is, at least in the outer parts of the core, it is composed of 
npe/^ matter. Hyperonic and quark matter may appear in the inner core without altering 
the results significantly as the binding energy is determined mainly by the physics of matter 
at saturation density and just above. We also assume that our selection of EoSs span the full 
range of EoS parameter space consistent with current observation. Regarding the last caveat, 
note that we are treating all other EoS properties (incompressibility, etc) as unconstrained 
- we have tested EoSs with as wide a range of those properties as possible. There are 
however, experimental constraints on some of those properties (e.g., the incompressibility of 
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symmetric nuclear matter at saturation density). Such constraints, and those that come in 
the future from experiment and observations, may improve the L—Mb relation and hence the 
constraints from L. Finally, we have taken the gravitational mass to be 1.2489M0 ignoring 
the observational error of ±0.0007Mq. This translates into an error of the same magnitude 
in Mb, i.e., of ~ 10^'^ Mq, which does not significantly change any of our conclusions. 



VI. CONCLUSION 



We have demonstrated that for a fixed gravitational mass Mq there exists a correlation 
between the binding energy, and hence the baryon mass Mb, of a neutron star and the 
slope of the nuclear symmetry energy sampled at around nuclear saturation density. The 
L — Mb correlation holds for a wide variety of EoSs calculated within a number of different 
theoretical frameworks and depends on only two assumptions: (1) that the star is indeed a 
neutron star and not a pure quark star (e.g. strange star) (2) that the EoS of the neutron 
star is sufficiently stiff at high densities to give a maximum mass above 1.44M0. Note that 
we do not exclude the appearance of hyperons, quark matter or other exotica in the denser 
parts of the star's core provided that they do not soften the EoS such that point (2) is 
violated. 

Combined with all current experimental data on L, this correlation gives an approximate 
constraint on the baryon mass of pulsar B of the double pulsar system PSR J0737-3039 
1.34 < Mb < 1.4M0. This is an independent constraint on the baryon mass to that derived 
from modeling the progenitor star to pulsar B under the assumption that it formed in an 
electron-capture supernova. We can expect the experimental constraints on L to become 
stricter in the next few years; the properties of isospin-asymmetric systems is at the forefront 
of experimental investigations at radioactive ion beam facilities (FAIR in Germany, RIKEN 
in Japan, SPIRAL2/GANIL in France, GSR in Ghina and FRIB in the United States) and 
the parity- violating electron scattering experiment at JLab (PREX), 36|. 

The e-capture constraints are shown to be consistent only with a value of L under ~70 
MeV, which falls in the lower half of the range constrained from a variety of nuclear ex- 
periments. If further experimental data suggests a value of L that is higher, the e-capture 
scenario may have to be revised, for example to include greater mass loss from the progenitor 
core during supernova, to shift the Mb constraint to lower values. The correlation explored 
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in this article reveals an intriguing new mutual interplay between the fields of nuclear ex- 
periment, nuclear theory, and astrophysical modeling, which we hope will encourage further 
experimental investigation and refinement of theoretical models in all areas. 
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FIG. 1: (Color online) The quantity RP ''•^^ versus neutron star radius K for the sample of EoSs, 
where the pressure P of neutron star matter is evaluated at saturation density no- 
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FIG. 2: (Color online) The pressure of /3-equilibrium matter at saturation density versus the slope 
of the symmetry energy L for the sample of EoSs. 
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FIG. 3: (Color online) The fractional gravitational binding energy of the neutron star BE/Mq 
versus the compactness parameter (3 for the sample of EoSs. The dashed line indicates the analytic 
Newtonian result. 
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FIG. 4: (Color online) The slope of the symmetry energy at saturation density L versus the baryon 
mass of a neutron star with a gravitational mass of 1.2489M0 for the sample of EoSs. The labora- 
tory constraints on L are represented by the shaded regions; the grey region 25MeV< L <115MeV 
is the most conservative constraint LI whereas the yellow region within is the constraint L2. The 
constraints on the baryon mass Mb from modeling pulsar B's progenitor under the assumption of 
an e-capture SN are given by the vertical lines; the dashed lines show the constraint Ml and the 
dash-dot lines the constraint M2. 
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